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Abstract 

The tradition in Navier-Stokes analysis of finding estimates in terms of the Grashof 
number Gr, whose character depends on the ratio of the forcing to the viscosity means 
that it is difficult to make comparisons with other results expressed in terms of Reynolds 
number Re, whose character depends on the fluid response to the forcing. The first task 
of this paper is to apply the approach of Doering and Foias \2'6\ to the two-dimensional 
Navier-Stokes equations on a periodic domain [0, L} 2 by estimating quantities of physical 
relevance, particularly long-time averages (•}, in terms of the Reynolds number Re = Ut/v, 
where U 2 = £~ 2 (||it||2) and t is the forcing scale. In particular, the Constantin-Foias- 
Temam upper bound pQ on the attractor dimension converts to a 2 Re (1 + lniZe) 1 ^ 3 , while 
the estimate for the inverse Kraichnan length is (a 2 Re) 1 / 2 , where ae is the aspect ratio of 
the forcing. Other inverse length scales, based on time averages, and associated with higher 
derivatives, are estimated in a similar manner. The second task is to address the issue of 
intermittency : it is shown how the time axis is broken up into very short intervals on which 
various quantities have lower bounds, larger than long time-averages, which are themselves 
interspersed by longer, more quiescent, intervals of time. 

1 Introduction 

1.1 General introduction 

In the last two decades the notion of global attractors in parabolic partial differential equations 
has become a well-established concept 0121131 H|- The general nature of the dynamics on the 
attractor A, in a time averaged sense, can roughly be captured by identifying sharp estimates 
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of the Lyapunov (or fractal or Hausdorff ) dimension of A, or the number of determining modes 
[3], with the number of degrees of freedom. Introduced by Landau [Bj, this latter idea says that 
in a dynamical system of spatial dimension d of scale L, the number of degrees of freedom TV is 
roughly defined to be that number of smallest eddies or features of scale A and volume X d that 
fit into the system volume L d 

A^(f) J . (1.1) 

This is the origin of the much-quoted M ~ Re 9 ^ result associated with the three-dimensional 
Navier-Stokes equations which rests on taking A ~ A& ~ LRe~ 3 ^, where A/% is the Kolmogorov 
length scale. In the absence of a proof of existence and uniqueness of solutions of the three- 
dimensional Navier-Stokes equations, at best this is no more than a rule of thumb result. It 
rests on a more solid and rigorous foundation, however, for the closely related three-dimensional 
LANS-a equations for which Foias, Holm and Titi [7j have proved existence and uniqueness of 
solutions. Following on from this, Gibbon and Holm [Sj have demonstrated that the dimension 
of the global attractor for this system has an upper bound proportional to Re 9 ^. An important 
milestone has been passed recently in another closely related problem with the establishment 
by Cao and Titi of an existence and uniqueness proof for Richardson's three-dimensional 
primitive equations for the atmosphere. 

For the Navier-Stokes equations the idea sits more naturally in studies in the two-dimensional 
context. The existence and uniqueness of solutions has been a closed problem for many decades 
and the nature of the global attractor has been well-established [1-5, 10-14]. While the two- 
and three-dimensional equations have the same velocity formulation, in reality, the former have 
a tenuous connection with the latter because of the absence of the drastic property of vortex 
stretching. As a result, the presence of vortex stretching in three dimensions, and perhaps 
other more subtle properties, have set up seemingly unsurmountable hurdles even on periodic 
boundary conditions. For problems on non-periodic boundaries, such as lid-driven flow, solving 
the two-dimensional Navier-Stokes equations is a technically more demanding problem - see 
some references in ^3 ^3 EES] ■ 

The sharp estimate found by Constantin, Foias & Temam 1 for the Lyapunov dimension of 
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the global attractor A expressed in terms of the Grashof number Gr 

d L (A) < ciGr 2 / 3 (l + lnO) 1/3 , (1.2) 

has been one of the most significant results in two-dimensional Navier-Stokes analysis on a 
periodic domain £2 = [0,L] 2 er . The traditional length scale in the two-dimensional Navier- 
Stokes equations is the Kraichnan length, which plays an equivalent role in two dimensions 
to that of the Kolmogorov length, A^, which is more important in three dimensions. In two 
dimensions, % and are defined respectively in terms of the enstrophy and energy dissipation 
rates e ens and e 



e ens = uL- 2 ^ \Vu\ 2 dVj , e = vLT 2 M \u\ 2 dVj , (1.3) 

where the pair of brackets (•) denote a long-time average defined as [2,3,10-13]. 

1 I* 

(#(")) = limlimsup- / g{r)dr . (1-4) 
9 (0) t Jo 

The inverse Kraichnan length rf^ 1 and the inverse Kolmogorov length A^~ 1 are defined in terms 



of e P „s and e as 



It has been shown by Constantin, Foias and Temam 1 that instead of using an estimate for 
e ens in terms of Gr, the upper bound for di can be re-expressed in terms of LrjZ 1 (see other 
literature on this topic [T71 HH1 

d L <c 2 {Lr)^) 2 {1 + In (L V ^)} 1/3 . (1.6) 

If dr, is identified with the number of degrees of freedom ]\[, this result is consistent with the idea 
expressed in (jl.lj) that in a two-dimensional domain, the average length scale of the smallest 
vortical feature A can be identified with the Kraichnan length rj^, to within log-corrections. The 
result in (|1.2|) has also been improved by Foias, Jolly, Manley and Rosa |2L)l I21j to an estimate 
proportional to Gr 1 / 2 (to within logarithmic corrections) provided Kraichnan's theory of fully 
developed turbulence is implemented |2*2] . 
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While these results display a pleasing convergence between rigorous estimates and scaling 
methods in the two-dimensional case, the tradition in Navier-Stokes analysis of finding estimates 
in terms of the Grashof number Gr, whose character depends on the ratio of the forcing to the 
viscosity v, means that it is difficult to compare with the results of scaling theories whose re- 
sults are expressed in terms of Reynolds number. One of the tasks of this paper is to estimate 
quantities of physical relevance, particularly long-time averages, in terms of the Reynolds num- 
ber, whose character depends on the fluid response to the forcing, and which is intrinsically a 
property of Navier-Stokes solutions. Doering and Foias [23] have addressed this problem and 
have shown that in the limit Gr — > oo, solutions of the d-dimensional Navier-Stokes equations 
must satisfy 1 

Gr <c(Re 2 + Re), (1.7) 

while the energy dissipation rate e has a lower bound proportional to Gr. The problem, however, 
is not as simple as replacing standard estimates in terms of Gr by Re 2 from Q1.7[) . Estimates 
such as that for di in Q1.2j) and the inverse Kraichnan and Kolmogorov lengths defined in 
ljl.51) . depend upon long time-averages of the enstrophy and energy dissipation rates defined in 
Other estimates of inverse length scales (to be discussed in ^1.2(1 also depend upon long 
time-averages. When estimated in terms of Re all these turn out to be better than straight 
substitution using Q1.7JI . These results are summarized in m.2l and worked out in detail in 

The second topic to be addressed in this paper is that of intermittency. Originally this 

important effect was considered to be a high Reynolds number phenomenon associated with 

three-dimensional Navier-Stokes flows. First discovered by Batchelor and Townsend |24| . it 

manifests itself in violent fluctuations of very short duration in the energy dissipation rate e. 

These violent fluctuations away from the average are interspersed by quieter, longer periods in 

the dynamics. This is a well established, experimentally observable phenomenon [25U261 l2"7j : its 

appearance in systems other than the Navier-Stokes equations has been discussed in an early and 

easily accessible paper by Frisch & Morf )28| . One symptom of its occurrence is the deviation 

of the 'flatness' of a velocity signal (the ratio of the 4th order moment to the square of the 2nd 
lr rhis result is not advertised in [33] but follows immediately from their equation (48). 
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order moment) from the value of 3 that holds for Gaussian statistics. 

Recent analysis discussing intermittency in three-dimensional Navier-Stokes flows shows that 
while it may be connected with loss of regularity, the two are subtly different issues j^H]- This 
is reinforced by the fact that although solutions of the two-dimensional Navier-Stokes equations 
remain regular for arbitrarily long times, nevertheless many of its solutions at high Re are known 
to be intermittent [313 [^ EH ESI EH • While three-dimensional analysis of the problem is based 
on the assumption that a solution exists |29l I35j . so that the higher norms can be differentiated, 
no such assumption is necessary in the two-dimensional case where existence and uniqueness 
are guaranteed. The result in both dimensions is such that the time-axis is broken up into good 
and bad intervals : on the latter there exist large lower bounds on certain quantities, necessarily 
resulting in their extreme narrowness and thus manifesting themselves as spikes in the data. 
This is summarized in <jl.2l and worked out in detail in 

1.2 Summary and interpretation of results 

For simplicity the forcing f(x) in the two-dimensional Navier-Stokes equations (divu = 0) 

u t + u • Vm = uAu - Vp + fix) (1.8) 

is taken to be divergence-free and smooth of narrow-band type, with a characteristic single 
length-scale t such that 1231135] 

I|v r 7ll2«r"||/|| 2 . (i.9) 

Moreover, the aspect ratio of the forcing length scale to the box scale is defined as 

a £ = L/e. (1.10) 

With frms = L~ d / 2 \\f\\2, the usual definition of the Grashof number Gr appearing in (|1.7|) in 
(^-dimensions is 




(1.11) 



The Reynolds number Re in Q1.7JI is defined as 

Re = — , U 2 = L- d {\\u\\l) , (1.12) 



19/08/06: arxiv.org/physics/0605086 Estimates for the 2D Navier-Stokes equations 6 

where (•) is the long-time average defined in (|1.4|) . One of the main results of this paper is the 
following theorem whose proof is given in M2.ll All generic constants are designated as c. 

Theorem 1.1. Let u(x, t) be a solution of the two-dimensional Navier-Stokes equations U.8\) 
on a periodic domain [0, L] 2 , and subject to smooth, divergence- free, narrow- band forcing f(x). 
Then estimates in terms of the Reynolds number Re and the aspect ratio an for the inverse 
Kraichnan length r/^ 1 , the attractor dimension di, and the inverse Kolmogorov length A^ 1 are 
given by 

Lr) k x <c(a\Refl 2 , (1.13) 
d L < cajRe [1 + In Re} 1/3 , (1.14) 

LX^ 1 <ca e Re 5/8 . (1.15) 

In the short proof of this theorem in ^12.11 the estimate for di in Q1.14JI is not re- worked from 
first principles but is derived from a combination of ()1.13j) and (j!.14f) . The result in (|1.15|) comes 
from a Re 5 / 2 bound on (Hi) and has also recently been found by Alexakis and Doering |36| . It 
implies that 

^ <ca t Re- x l 2 , (1.16) 

whereas in three-dimensions the right hand side is 0(1). The estimate in (|1.14|) is also consistent 
with the result of Foias, Jolly, Manley and Rosa |20| when their Gr 1 ^ 2 estimate is converted to 
one proportioanl to Re. Their estimate, however, was based on the implementation of certain 
features of the Kraichnan model |22| . while (|1.14|) is true for all solutions and requires no 
assumption of fully developed turbulence. 

The estimates for r/^ 1 and dr, are consistent with the long-standing belief that Re 1 ! 2 x Re 1 ^ 2 
grid points are needed to numerically resolve a flow; indeed, when the aspect ratio is taken into 
account, Theorem II. II is consistent with aiRe 1 ^ 2 x aiRe 1 ^ 2 . However, both these estimates are 
dependent upon only the time average of low moments of the velocity field. For non-Gaussian 
flows, low-order moments are not sufficient to uniquely determine the statistics of a flow. Thus it 
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is necessary to find ways of estimating small length scales associated with higher-order moments. 
In ^2, 21 we follow the way of defining inverse length scales associated with derivatives higher than 
two, introduced elsewhere ^H], by combining the forcing with higher derivatives of the velocity 
field such that 

F n = f (|v n w| 2 + r 2 |V n /| 2 ) dV, (1.17) 
Jn 

where \Gr{\ + In Gr)] 2 is a characteristic time: this choice of r is discussed in 

Appendix The gradient symbol V n within (|1.17|) refers to all derivatives of every component 
of u of order n in L 2 (Q). The F n are used to define a set of time-dependent inverse length scales 



(77 \ 2(n-r) 
TJ ■ (L18) 

Actually, k 2 ™^ behaves as the 2nth-moment of the energy spectrum as shown by 

More relevant to the two-dimensional case, fO™ ^ behaves as the 2{n — l)th-moment of the 
enstrophy spectrum. Using Landau's argument the dimension of the global attractor di(A) 
was identified with the number of degrees of freedom J\f. In a definition was introduced to 
represent the number of degrees of freedom associated with all higher derivatives of the velocity 
field represented by /i„ jr , which is itself an inverse length. This naturally leads to the definition 
of the infinite set 

Mn, r = L 2 ( K l r ) . (1.20) 

Using the definition of the quantities A nj o an d A n> i (n > 2) 

3n-2 3n-4 , 

A„o = , A„i = -, (1-21) 

2n ' ' 2(n - 1) K J 

the second main result of the paper is a theorem whose proof is given in H2.2I : 

Theorem 1.2. Let K n>r be the moments of a two-dimensional Navier-Stokes velocity field defined 
in U.18\) . Then in a two-dimensional periodic box of side L the numbers of degrees of freedom 
A/" n ,i and M n fi defined in il.20\) are estimated as (n > 2) 

M n ,i < Cn^ajRe)^- 1 (1 + In Re) 1 ' 2 , (1.22) 
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Mn,0 < Cn,o(a| 



1/2 



(1.23) 



where A ni o and A n) i are defined in M.21\) . 

Note that A2 i o = A2 i = 1. Thus the estimate for the first in each sequence, A/2,1 and M\ Oi are of 



logarithm. The exponents in (|1.22j) and (|1.23l) provide an estimate of the extra resolution that 
is needed to take account of energy at sub-Kraichnan scales. Notice that in the limit n — > 00 
both exponents converge to 3/2. 

The intermittency results of show that there can exist small intervals of time where there 



of spiky data whose duration must be very short. Estimates are found for the width of these 
spikes which turn out to be in terms of a negative exponent of Re. 

2 Time average estimates in terms of Re 
2.1 Proof of Theorem fTTTl 

The first step in the proof of Theorem ll.il which has been expressed in HI. 21 is to find an upper 
bound on (H2) in terms of Re. Consider the equation for the two-dimensional Navier-Stokes 
vorticity u = lo k 



the same order as the estimate for dx, namely ojRe(l + In Re) 1 / 3 except in the exponent of the 





(2.1) 



and let H n be defined by (n > 0) 




(2.2) 



For a periodic, divergence-free velocity field u 




(2.3) 
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Then the evolution equation for Hi is 

iHi = -vH 2 + [ u-cmlfdV (2.4) 

< -vH 2 + \\u\\ 2 \\V 2 f\\ 2 (2.5) 

< - I /i? 2 +r 2 ||w|| 2 ||/|| 2 , (2.6) 

where the forcing term has been integrated by parts in 1)2. 4 j) and the narrow-band property has 
been used to move from Q2.5|) to (|2.6|) . Using the definitions of Re, Gr, and ag in (|1.12l) . (|1 . Ill) 
and (|1.1U|) , the long-time average of H 2 is estimated as 

(H 2 ) < L 2 r Q v 2 ReGr (2.7) 

< cajr 4 u 2 Re 3 + 0(Re 2 ). (2.8) 

This holds the key to the three results in Theorem ll.il 

The inverse Kraichnan length — e ens jv 3 with e ens = vL 2 (H 2 ), can now be estimated 
by noting that 

L 6 e ens < ca\v 3 Re 3 (2.9) 

and so 

Lr]- 1 <c(a 2 e Re) 1/2 , (2.10) 

which is (|1.13|) of Theorem (|1.1|) . The estimate for d,L in (|1.14|1 then follows immediately from 
the relation between the estimate for in (|1.6|) and (|2.1U|) . 

Finally, we turn to proving the estimate for in (|1.15|) which turns around the use of the 
simple inequality H 2 < H 2 Hq. The next step is to use the fact that 

{Hi) < {H 2 f 2 (H ) 1/2 (2.11) 

= va e Re (H 2 } 1/2 . (2.12) 

Using the upper bound in (|2.7() gives 

(Hi) < cu 2 a 2 e e 2 Re 5/2 , (2.13) 

which then gives (jl.l5|) in Theorem (Jl.lj) . In fact, ()2.13|) is an improvement in the bound for 
(Hi) from Re 3 to i?e 5//2 . This result has also been found recently by Alexakis and Doering 36 . 
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2.2 Proof of Theorem O 

Having introduced the notation for H n in (|2.2[) . similar quantities are used that contain the 
forcing |351 I29| . namely 

F n = I (|V n u| 2 + r 2 |V n /| 2 )(iy, (2.14) 
Jn 

defined first in (|1.17|) . and the moments K n ^ r defined in (|1.18() 

/ P \ 2(n-r) 

n n , r (t) := . (2.15) 



The parameter r in (|2.14|) is a time scale and needs to be chosen appropriately. The idea is 
that it should be chosen in such a way that the forcing does not dominate the behavior of the 
moments of the velocity field. Defining ujq = £~ 2 v, it is shown in Appendix [S] that this end is 
achieved if t~ 1 is chosen as 

r- 1 = uj [Gr(l + In Gr)} 1 ' 2 (2.16) 
< cuJ Re(l +lni?e) 1 / 2 . (2.17) 

As a preliminary to the proof of Theorem 11.21 we state the ladder theorem proved in |35l I29j . 
Theorem 2.1. The F n satisfy the differential inequalities 

i Fo < - v f 1 + ct- 1 Fq, (2.18) 

\F X < -vF 2 + ct~ 1 F 1 , (2.19) 

and, for n > 2, either 

\F n < -z.F n+1 +c ni i(||Vit|| 00 + T- 1 )F n , (2.20) 

\F n < -lvF n+1 + c^fV 1 !!^^ + T-^F n . (2.21) 

The L°°-inequalities in Theorem 12.11 particularly HVwjloo in (|2.20|) . can be handled using a 
modified form of the L°°-inequality of Brezis and Gallouet that has already been proved in ^Hj : 



or 
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Lemma 2.1. In terms of the F n of {2.1$ and Ksp °f i<i.l5\) . a modified form of the two- 
dimensional L°° -inequality of Brezis and Gallouet is 

IIVuHoo < cFl 12 [1 + ln(L K3 . 2 )] 1/2 . (2.22) 

This lemma directly leads to an estimate for (k 2 r ) for r > 2. 
Lemma 2.2. For n > r > 2, to leading order in Re, 

L 2 (4, r ) < c {a 2 e Ref /2 {l + In Re) 1/2 . (2.23) 
Proof: By dividing (|2.2Uj) by F n and time averaging, we have 

v(i&+i, n ) < c n ,i(||Vw|| 00 ) + co;o^e(l + ln J Re) 1 / 2 . (2.24) 

However, because K nx < K n+ i >n for r < n, for every 2 < r < n, in combination with Lemma 12. 11 
we have 

v <K 2 ir ) < c (V 2 1/2 [1 + ln(LK 3j2 )] 1/2 ) + cuj Re{l + In Re) 1 / 2 . (2.25) 
The logarithm is a concave function and k^2 < K n ,r so Jensen's inequality gives 

L 2 « r ) < L 2 i/ -1 c (F 2 ) 1/2 < [1 + ln{L 2 « r }}] > 1/2 + ca 2 Re{\ + In fe) 1 ^ . ( 2 .26) 

The estimate for (F2} can be found from {H2) in (|2.7|1 : the extra term t 2 ||V 2 / U 2 . is no more than 
0{Re 2 ). Standard properties of the logarithm turn inequality l)2.26p into ()2.23j) . ■ 

Lemma 12.21 gives estimates for (k 2 r ) for r > 2. These are used in the following theorem to 
give better estimates for the cases r = and r = 1. Prior to this, it is necessary to state the 
results that immediately derive from Q2.18JI and (|2.19l) by respectively dividing through by Fq 
and F\ before time averaging 

M,o = l2 ( k i,o) < ca}Re{l + \nRe) 1/2 , A/2,1 = L 2 < ca 2 e Re(l + lnRe) 1/2 . (2.27) 

With the estimates in (|2.27|) we are now ready to complete the proof of Theorem 11.21 
Proof of Theorem 11.21 Let us return to (|2.23|) in Lemma l2~2*l and use the fact that 
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and thus 

n — 2 1 

Kl)^^)"- 1 ^,!)"" 1 • ( 2 - 29 ) 
Using (|2~2~3|) in Lemma E21 together with (f2~2Tj) . for n > 2, 

Mn,i = L 2 (kI a ) < c n ,i (ajRe)^^ [1 + lni?e] 1/2 , (2.30) 

which coincides with a 2 Re(l + ln Re) 1 / 2 at n = 2 but converges to i?e 3//2 (l + ln Re) 1 / 2 as n — > oo. 
The exponent A rai i is defined in (|1.21|) . 

Likewise, in the same manner as (|2.28j) we have 

«o><«i> n Ko)" ■ (2-31) 

Thus we find that for n > 2 

K,0 = L 2 « > < c nfi (a}Re)^r [1 + lnRe] 1/2 . (2.32) 

The exponent A n is defined in (jl.21|) . ■ 

3 Point-wise Estimates 

Let us consider the differential inequalities for Hq and H\: 

^K-vHr + WfhH] 12 , (3.1) 

\Hx < -uH 2 + r 2 ||/|| 2 F 1/2 , (3.2) 

having used the narrow-band property on (|3.2j) . Upon combining Poincare's inequality with 
Lemmas IB . 1 1 and IB . 21 in Appendix [B] we obtain 

amt^ooHo <ca\v 2 Gr 2 < c a\v 2 Re A , (3.3) 

and 

lim^oofli < cr 2 4v 2 Gr 2 < C r 2 aju 2 Re 4 . (3.4) 

The additive forcing terms in F\ and Fq are of a lower order in Re so we end up with 

lim^ooFo < ca 6 e v 2 Re 4 + 0(Re 2 ) , (3.5) 
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lim^oo-Fi < ct" 2 a\v 2 Re 4 + 0{Re 2 ) . (3.6) 

The estimate for F\ enables us to obtain point-wise estimates on F n , n > 2 18; ( sec. 7.2]. In 
fact we have the following lemma. 

Lemma 3.1. As Gr — > oo 

lim^ooi 7 ; < c n v 2 r 2n al n Re 4n . (3.7) 
Proof: Applying a Gagliardo-Nirenberg inequality in two-dimensions to Vw we obtain 

II Vu||oo < c || V n it||^||Vu||^ a < cElF^ , (3.8) 
with a = Using this in (|2.20f) gives 

\F n < -vF n+l + c n Fn + ^F^ + cuj Re(l + \nRef' 2 F n . (3.9) 
Moreover the following inequality can easily be proved using Fourier transforms 

F P N q < n~ p F p N+q , (3-10) 
from which, with N = n, p = n — 1, q = 1, it can be deduced that 

n 

pn-l 

-F n+ i < V . (3.11) 

F n-1 

We now use ()3.11|) in ()3.9|) to obtain 

n 

F"- 1 



F n <-v^^ + c n F n +2 F l 2 +cuj Q Re{l + \nRe) 1 / 2 F n , (3.12) 



i 

2-l n 



fn-1 



with a = zrzT- We use now estimate (j3.6j) in (|3,12j) with the further use of Lemma lB.2l to obtain 

ImT^oo^n < c n v 2 t- 2n afGr 2n , (3.13) 
which leads to the result. ■ 

The above Lemma enables us to obtain an estimate on the wave-numbers K nr . 
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Lemma 3.2. For n > r > 0, as Gr — > oo 



lim t _ ) . 0O (LK n)r ) < c n a £ " r Re n ~ r (1 + lnRe) 2> - n - 



(3.14) 



Proof: Essentially one uses the upper bound on F n and the lower bound on F r which can be 
calculated from the forcing part in terms of Gr, leading to the result (see also [181 Ch. 7]). ■ 

4 Intermittency: good and bad intervals 

The issue of intermittency in solutions of the two-dimensional Navier-Stokes equations is now 
addressed. While the F n and bounded from above for all time, nevertheless it is possible 

that their behaviour could be spiky in an erratic manner. To show how this might come about, 
consider the definition of K n ^ T in (|1.18|) from which we find 



P _ Ji ( K n+l,t 



2(n+l-r) 



F, 



(4.1) 



Now consider inequality (|3.9j) re-written as 

2rc 



< —UK 



n,l 



+ On 



Kn+1,1 



1/2 

Kn,iF 1 + c loq Re (I + In Re) 



1/2 



(4.2) 



Fn. V ^n,l / V 

where we have used (|4.1|) and the fact that n n i < n n +i,i m the middle term. Using Young's 
inequality on this same term we end up with 

p , \ 2n 



i - n ^ i 2 / K n+l,l 



^±hl) +c n v- 1 F 1 +cuj Re(l + In Re) 1 ' 2 . (4.3) 

«n,l / 

The main question is whether, for Navier-Stokes solutions, the lower bound on 



> 1 



(4.4) 



can be raised from unity. A variation on the interval theorem proved in is used. 



Theorem 4.1. For any value of the parameter fx 6 (0, 1), the ratio K n +l i/ K ni obeys the long- 
time averaged inequality (n > 2) 



(i 2 <l)^ 



{a 2 e Re) A ^{l + mRey/ 2 



> 0. 



(4.5) 
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where the c n are the same as those in Theorem M.fyl . Hence there exists at least one interval of 
time, designated as a 'good interval', on which the inequality 

^ 2 o.n a. 77 , ,_ „.m/9 ( 4 ' 6 ) 



«n,l / ~ {a}Re) h ^(l+\nRe) l l 2 
holds. Those other parts of the time-axis on which the reverse inequality 

C„ — < ^ ; 



V Kn J (ajRe) A ^(l + liiRe) 1 / 2 K ' ' 

holds are designated as 'bad intervals'. 

Remark: In principle, the whole time-axis could be a good interval, whereas the positive time 
average in Q4.5J) ensures that the complete time-axis cannot be 'bad'. This paper is based on the 
worst-case supposition that bad intervals exist, that they could be multiple in number, and that 
the good and the bad are interspersed. The precise distribution and occurrence of the good/bad 
intervals and how they depend on n remains an open question. The contrast between the two- 
dimensional and three-dimensional Navier-Stokes equations is prominent; while no singularities 
can occur in the k Uj i in the two-dimensional case, in three dimensions it is within these bad 
intervals that they can potentially occur. 

Proof: Take two parameters < \i < 1 and < ct < 1 such that fj, + a = 1. The inverses fi^ 1 
and a -1 will be used as exponents in the Holder inequality on the far right hand side of 

2a \ I / \ 2a/ n\ P 

K n+l,l \ 2a \ ^ I I K n+l,l \ \ / 2 



thereby giving 

(4.9) 




/ / K n+lt 
\ \ K n,l 

Two-dimensional Navier-Stokes information can be injected into these formal manipulations: 
the upper bound on j) from Theorem (|1.2|) and the lower bound Ln n) \ > 1 are used in the 
ratio on the far right hand side of ()4.9j) to give (|4.5j) . with the same c n as in Theorem (jl.2|) . ■ 
Now consider what must happen on bad intervals. It is always true that «n+l,i/ re ni > 1, so 
(|4.7j) implies that on these intervals there is a lower bound 

L\ 2 n 1 > c ri {a}Re) k ^^{l + \nRe) 1 ^ . (4.10) 
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This lower bound cannot be greater than the upper point-wise bound in (|3.14|) . which means 
that // is restricted by 

An, < 2 (2n-l , _ (4n) 



\i \ n — 1 

Moreover, the factor of in the exponent makes the lower bound in (|4.1Uj) much larger than 
the upper bound on the average (k^ x ) given in Theorem H1.2|) . These intervals must therefore 
be very short. To estimate how large they can be requires an integration of Q4.3J) over short 
times At = t — to which, in turn, requires the time-integral of H\ for short times At. We use 
the notation f^ t = , with the definition ojq = vl~ 2 . 



Lemma 4.1. To leading order in Re 

I F\dt < va\ [ciaj + C2 woAt] Re 4 . 
J At 

Proof: Integrating (|3,lj) over a short time At gives 

v I H x dt < ±H (t ) + At[£- 2 v 3 afGr 2 ] 
J At 

< cxa$v 2 Re 4 + At [c 2 £' 2 u 3 ajRe 4 ] , 
having used (|3,3|) for the i-£fo(^o)-t erm - The forcing term in F\ is only 0(Re* 



(4.12) 



(4.13) 



«n,l(t) 
max(K ni i 

















/./. ,. i > Re^/f ! 






I 


\ i 
\ i 

\ 1 




■ 






Lo^^l^^ 







A 



(At) b (At)g 

Figure 1: A cartoon, not to scale, of good/bad intervals for some value of n > 3. 
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Now we wish to estimate uj^At in terms of Re. Integrating (|4.3|) . using (|4.13|) and the lower 
bound (|4.1Uj) and multiplying by £ 2 , we have 

\i 2 [In F n (t)- In F n (to)] + \c n va- 2 {ajRe)^^{l + In Ref^At 

< t 2 a\ [cia| + c 2 uj Q At] Re A 

+ c£ 2 u AtRe(l + \nRe) 1/2 . (4.14) 

As Gr — > oo, the dominant terms are 

w A£ {a7 2 (a|i?e) A - 1 ^(l + In Re) 1 ^ - a^e 4 } < Cl a\Re A . (4.15) 
Choosing ji in the range, to leading order we have 

H < \K,i , (4.16) 

then At must satisfy 

uj At < c{ajRef- K ^^. (4.17) 

Because the exponent in (|4.17j) is necessarily negative these intervals are very small and decreas- 
ing with increasing Re. Combining Q4.11J1 with (|4.16|) we have 

(n — 1) 

5 L_2 y A B , 1 < / x<iA B4 , (4.18) 

which actually holds for every n > 1. Figure 1 is a cartoon- like figure displaying the lower 
bound on the bad intervals of width (At)b and also the maximum of K n i allowed by (|3.14|) in 
Lemma 13.21 The full dynamics of two-dimensional Navier-Stokes is actually determined by the 
intersection of all cartoons for every n > 3 on the grounds that the position and occurrence of 
the bad intervals varies with n. Thus we are interested in the limit n — > oo which determines 
that the range of [i is squeezed between 

l(>-l)<»<l(i-l)- 

Thus, in the limit, fi takes a value just under 3/8. We conclude that the interval theorem 
(Theorem 14. 1|) reproduces the effects of intermittency in a two-dimensional flow by manifesting 
very large lower bounds within bad intervals and suppressing spiky behaviour within the good 
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intervals which must be quiescent for long intervals, otherwise the long-time average would be 
violated. 
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A Forcing & the fluid response 

For technical reasons, we must address the possibility that in their evolution the quantities H n 
might take small values. Thus we need to circumvent problems that may arise when dividing by 
these (squared) semi-norms. We follow Doering and Gibbon 35j who introduced the modified 
quantities 

F n = tf n + r 2 ||V n /|||, (A.l) 

where the "time-scale" r is to be chosen for our convenience. So long as t / 0, the F n are 
bounded away from zero by the explicit value t 2 L 3 £ n f 2 ms . Moreover, we may choose r to 
depend on the parameters of the problem such that {F n ) ~ (H n ) as Gr — > oo. To see how to 
achieve this, let us define 

t = £ 2 u- 1 [Gr(l + In Gr)}- 1/2 . (A.2) 
Then the additional term in (|A.1|) is 

r 2 ||Wlll = tfv-H^f^Gril + lnGr)]- 1 

= v 2 r {2n+2 ^U i Gr{l + \nGr)- 1 . (A.3) 

Now Doering & Foias [23] proved that in d-dimensions, the energy dissipation rate e has a lower 
bound of the form 

e>cv 3 r z L' x Gr. (A.4) 
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Using this on the far right hand side of (jA.3j) we arrive at 

r 2 ||V n /||| < c 6 e r( 2n -V L^is- 1 (1 + In Gr)- 1 
f L s (2n-l) 



j J L-^-V^il + mGr)- 1 . (A.5) 



Using Poincare's inequality in the form Hi < (27rL) 2 ( n ^H n , as Gr -^oowe have 

,2||y7nf||2 

' ^{" 2 < c 6 af "- 1} (1 + In Gr)- 1 . (A.6) 

\ H n) 

Hence, the additional forcing term in (|A.1|) becomes negligible with respect to (H n ) as Gr — > oo, 
so the forcing does not dominate the response. 

B Comparison theorems for ODEs 

We present a comparison theorem for ODE which is useful for obtaining various estimates. We 
start with the following classical result. 

Lemma B.l. Let f : [0, T] xR — > R be a continuous function which is locally Lipschitz uniformly 
in t: for all intervals [a,b] C R there exists a constant such that \ f(s,x) — f(s,y)\ < C\x — y\ 
for all x, y € [a, b] and all s € [0, T]. Furthermore, let x G j4C([0, T],R) be such that 

x(t) < f(t,x(t)) 

for all t E [0, T] and let y(t) be the solution of y(t) = f(t,y(t)) on [0,T]. Assume further that 
x(0) < y(0). Then, x(t) < y{t) for all t G [0,T]. 

We can use this Lemma to prove the following useful result. 

Lemma B.2. Let x : [0, T] — > [0, oo) be an absolutely continuous function with x(0) > which 
satisfies 

x < A x + Fx ni -Ex n \ (B.l) 
where Ao, F, E > and 1 < n\ < rii. Then 

limsupx(t) < (AAqE- 1 )^ 1 + (2FE~ 1 )^^. (B.2) 
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